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Effective spin Hall properties of a mixture of materials with and without spin-orbit 
coupling: Tailoring the effective spin-diffusion length 
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We study theoretically the effective spin Hall properties of a composite consisting of two materials 
with and without spin-orbit (SO) coupling. In particular, we assume that SO material represents 
a system of grains in a matrix with no SO. We calculate the effective spin Hall angle and the 
effective spin diffusion length of the mixture. Our main qualitative finding is that, when the bare 
spin diffusion length is much smaller than the radius of the grain, the effeetive spin diffusion length 
is strongly enhanced, well beyond the “geometrical” factor. The physical origin of this additional 
enhancement is that, with small diffusion length, the spin current mostly flows around the grain 
without suffering much loss. We also demonstrate that the voltage, created by a spin current, is 
sensitive to a very weak magnetic field directed along the spin current, and even reverses sign in a 
certain domain of fields. The origin of this sensitivity is that the spin precession, caused by magnetic 
field, takes place outside the grains where SO is absent. 

PACS numbers: 85.75.-d,72.25.Rb, 78.47.-p 


I. INTRODUCTION 

The spin Hall effecfP^S HE), predicted theoretically 
more than four decades agcP^, is nowadays routinely 
observed in many materials,®^ which include tradi¬ 
tional and exotic metals, prominent semiconductors, and 
graphene. Moreover, the inverse spin Hall effect (ISHE), 
i.e. generation of voltage drop normal to the spin current, 
was recently “put to work”. It serves as a tool to detect 
whether or not the spin current is injected into a non¬ 
magnetic material from an ac-driven ferromagnet in the 
course of spin pumping. Most recentl^ff^H^ the pumped 
spin currents in certain polymers were registered via in¬ 
verse spin Hall voltage which they induced in Pt elec¬ 
trode located at some distance from the interface with 
ferromagnet. 

The latest focu^^^^^of the research on the spin physics 
in organics is the study of the properties of platinum- 
containing TT-conjugated polymers. In these materials Pt 
atoms are embedded in the polymer backbone chains. 
While the SO coupling, which is the origin of the SHE, is 
very weak in polymers, adding of Pt creates the elements 
of the backbone where it is locally strong. These elements 
can be separated either by one or by three 7r-conjugated 
spacer unit lengths. In this regard, a general question 
arises: how the spin Hall effect is realized in composite 
materials where the strong SO and low SO domains are 
intermixed? Note that, by now, all theoretical studies 
of SO-related transport assumed that the SO coupling is 
homogeneous. 

The goal of the present paper is to develop an ele¬ 
mentary theory which addresses the question formulated 
above. Unlike Refs. [25] and [26l we will not specify a 
mechanism of SO on the microscopic level, but rather 
focus on purely “geometrical” aspects. Namely, we will 
consider the following minimal model: a system of SO 
grains is dissolved in a matrix with no SO. The question 



FIG. I. (Color online) (a) Conventional geometry for the 
inverse spin Hall effect. Spin current flowing along y causes 
a buildup of the voltage, between the edges z = -EL 12. 

The buildup takes place as long as y is smaller than the spin 
diffusion length, A. (b) Schematic illustration of a “granular” 
geometry, where the SO-coupled material is dissolved in the 
matrix with no SO coupling, (c) microscopic scenario of ISHE 
on a single spherical granule of a radius, a. Spin current with 
polarization along x turns the sphere into an electrical dipole 
directed normally to the current. The magnitude of a dipole 
moment, Vc, depends on the ratio between a and A, while the 
electric field inside the granule is homogeneous. 


we will be interested in is: what are the effective spin 
Hall characteristics of the mixture. 

Eirstly, we address a mechanism of the formation of 
the inverse spin Hall voltage between the edges of the 
sample in the geometry of the mixture. Unlike the case 
of homogeneous SO, this formation happens as follows. 
The spin current turns each SO grain into an electric 
dipole. All dipole moments are oriented normal to the 
spin current. Thus the potentials they create at the up¬ 
per and the lower boundaries of the sample add up. The 
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difference of these potentials is the effective ISH voltage, 
Kff 5 of the mixture, which can be related to the effective 
spin Hall angle, 0 ^^. 

Naively, one would expect that, in a mixture of grains 
of density, n, and radius, a, the relation 0 ^^ = {na^) 0 ^^ 
holds within a numerical factor. Here 0 ^^ is the spin Hall 
angle of the bulk SO material. This is simply because 
na^ is the volume fraction of the SO material. Equally, 
one would expect that the effective spin relaxation time 
of the mixture is l/(na^) times longer than in the SO 
material, so that spin diffusion length, Aeff is related to 
the spin diffusion length. A, of the SO material as Aeff = 

The above expectations are correct only in the limit 
when the grains are small enough, namely, a ^ A, so 
that the portion of spin polarization, which is lost within 
a single grain, is small. The opposite case of large grains, 
a ^ A, is much less trivial. As we show below, in this 
limit ^ Ana^E®^, while Aeff ^ other 

words, at small A, the effective spin-diffusion length sat¬ 
urates. This finding can be loosely interpreted from the 
perspective of diffusion in the presence of the absorbing 
traps. The stronger is the absorption, the smaller is the 
concentration of particles at the position of the trap. 

Finally, we will demonstrate that EJf is sensitive to a 
very weak magnetic field. In a homogeneous material, the 
spin Hall effect gets suppressed in the field with Larmour 
frequency ft ^ where is the spin-relaxation time. 
For the mixture, the characteristic field is ^ T“^, where 
T is the diffusion time between the sample edges. This 
is because spin precession takes place mostly outside the 
grains. The paper is organized as follows. In Sect. H 
we solve an auxiliary problem of electric the polarization 
of a grain with a given radius, a, by the spin current. 
The solution is then employed to calculate the effective 
inverse spin Hall voltage in the mixture of grains with 
concentration, n. Sensitivity of this voltage to a weak 
longitudinal magnetic field is studied in Sect. HI. In Sect. 
IV the effective diffusion length, Aeff, of the mixture is 
expressed via A, a, and the parameter na^. The physics of 
elongation of Aeff for small A ^ a is discussed in Sect. V. 
Concluding remarks are presented in Sect. VI. 


II. CALCULATION OF EFFECTIVE 
CHARACTERISTICS OF THE MIXTURE 

A. Single grain 


The simplest way to incorporate the spin Hall effect 
on a quantitative leveP^ is to add to the current density, 
j = aE, the term ^Dcm\P where a and D are the 
conductivity and the diffusion coefficient, respectively, 
P{r) is the coordinate-dependent spin polarization. The 
strength of the SO coupling is quantified by a dimension¬ 
less parameter 7. The system of coupled equations for 



FIG. 2. (Color online) (a) The cross section z = 0. Distri¬ 
bution of the spin current in the {x, y)-plane in the presence 
of a spherical grain, Eq. (38), is illustrated schematically for 
r > a. Inside the grain, r < a, this distribution is determined 
by Eq. ( [2^ . (b) Distribution of the spin polarization along 
the radius, r, is plotted for Dout/Din = 2 and three values of 
A: a/\ — 0.2 (green), a/\ — A (blue), and a/\— 12 (purple). 
Enhancement of the effective spin diffusion length for small 
A/a is a result of a strong suppression of polarization near the 
boundary r — a. 


the spatial distribution of P{r) and j{r) read^ 


j = cfE e^DcnrlP. (1) 

dP V 

Qij = —o’SijkPk‘ ( 2 ) 

The second equation defines the component i of the flux 
of the .7-pro jec tion of spin polarization. The system be¬ 
comes closed^ when it is complemented by the continuity 
equation 


dQij Pj _ 
dx, T, 


(3) 


Consider an isolated spherical grain with radius, a, and 
with the strength of S0-coupling, 7, embedded into an 
infinite medium with 7 = 0 and with no spin relaxation, 
r, = oc. Fig. B Assume that the flux of spins, oriented 
along the x-axis and flowing along the 7/-axis, is incident 
on the grain. In application to the geometry. Fig. the 
essence of the inverse spin Hall effect is that the incident 
spin current, induces an effective electric dipole on the 
sphere. The induced dipole moment is perpendicular to 
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both, the current direction and polarization direction in 
the incident flux, i.e. it is directed along the 2 ;-axis. 

To calculate the magnitude, Vc^ of the dipole moment 
it is natural to switch to spherical coordinates in which 
the incident polarization, and the spin- 

current density, iy = is ^ have the form 


'i 

P = — ^ r sin 6 ^e 0 , (sin + cos 6 ^ 60 ), (4) 

-^out 

where e^, e^, and are the unit vectors along radial, 
polar, and azimuthal axes, respectively, see Fig. 

Induced dipole moment along 2 : creates an electrostatic 
potential. 


cos 0 

out ^2 5 


(5) 


outside the sphere. 

From the form of we conclude that the 0- 

dependence of cp inside the sphere is also proportional 
to cos 6 >. This, together with Poisson’s equation Aip = 0, 
suggests that the induced electric field, inside the 
sphere is homogeneous, so that 


<fin =-Ei^r cose. ( 6 ) 

Substituting Eq. ([^ into Eq.(|^, and taking into ac¬ 
count that dE^^jdXi = 0 , we conclude that all the com¬ 
ponents of polarization inside the sphere satisfy the dif¬ 
fusion equation 

A„AP,- + ^ = 0, (7) 

where is the diffusion coefficient inside the sphere. 

As we will see below, the polarization, P{r), has only 
^-component inside the sphere and at all distances out¬ 
side the sphere. As in the incident flux, Eq. Q, the 
angular dependence of is oc sin^. Outside the sphere, 
where AP = 0, the general form of P^ is 

Pout = sin 6 >e^, ( 8 ) 

where the constant is the “spin polarizability”. Inside 
the sphere, the solution of Eq. proportional to sin^, 
has the form 


Pi„ = P ii(r/A)sin0e^, (9) 


where P is a constant, and 

A = , (10) 


is the diffusion length. The function ii{x) is a modified 
spherical Bessel function. We chose the function ii be¬ 
cause it is finite at x = 0 . 

While the polarization has only 0-component, the spin 
current, defined as a flow of the ^-component of spin, can 
be presented in the vector form 




is (1 


2 ^ 

J.3 


) sin^e^ + (1 + 


Xs \ o 

^jcost/e^J, 


( 11 ) 


where the first term is dP^^^/dr^ while the second term 
is (^) dP^^JdO. At large distances the current Eq. ( 11 ) 
reproduces Eq. 

There are two unknown constants, P^ and y^, in the 
expressions for electric field and spin polarization inside 
the sphere, and two unknown constants, and P, in 
the corresponding expressions outside the sphere. These 
constants are determined from the four boundary condi¬ 
tions at r = a: 

(i) Continuity of the tangent component of electric field 




3 ■ 


( 12 ) 


(ii) Continuity of the normal component of the charge 
current 


267A„P. , 2(T„„tP, 

H- ii{a/\) = -^— 


(13) 


(iii) Continuity of the spin polarization 


(a+^) =Pii(a/A). (14) 


Po 

(iiii) Continuity of the spin flux though the boundary 


DP ^ 


(t - ')■ O'i) 


The system Eqs. -(p!^ yields the sought expression 
for the spin-current-induced dipole moment 






((7i„ + 2(7o„t)A4 


(16) 


where A4 in denominator is the dimensionless combina¬ 
tion 


M = 


2D„ 


27^0-in 


A, 


( 7 i„ -h 2 c 7 o 


Ai(f) 


(17) 


Naturally, the proportionality coefficient between P^ and 
the spin current contains the first power of the SO cou¬ 
pling strength, 7 . 

The second term in Eq. ( [T7| ) contains 7 ^, and can 
be safely neglected. The ratio can be replaced 

by cTout/c^in- It is seen from Eq.(pT|) that the factor Ad 
depends strongly on the relation between the radius of 
the sphere and the spin-diffusion length. Eor a X the 
last term in Eq. ( pT[ ) is 1, while for A ^ a it is big and 
equal to a/A. In the latter case Eq. (16) yields P^ oc Xa^. 
This dependence has a simple interpretation. Namely, 
for A ^ a the induced dipole is generated only inside a 
spherical layer of a thickness ^ A near the surface of the 
sphere, see Eig. 

Description of a direct spin Hall effect for a sphere is 
completely similar to the case of the inverse spin Hall 
effect considered above. A charge current, ic, along the 
y direction generates a spin dipole moment, P^, in the 
^(-direction. Analytical expression for P^ is similar to Eq. 
( 116 ) 


Vs = 


3(Ti„a^7 


e(cri„ -h 2(7out) A„A4 


(18) 




















4 



FIG. 3. (Color online) The effective spin Hall voltage is 
the sum of contributions from individual SO-induced dipoles. 
With density of granules, n, the typical distance between the 
neighbors is It is much bigger than the radius, a, but 

much smaller than the sample width, L, which allows to re¬ 
place the sum by the integral Eq. (20). 


FIG. 4. (Color online) Dependence of the effective ISHE volt¬ 
age on a longitudinal magnetic held, is plotted from Eq. 
(25) for three different positions, yo, in the units of (dL)^/^, 
along the sample. Blue, violet, and green curves correspond 
to yo = 2, yo =4, and yo = 6 , respectively. 


reduces to multiplication by d. The final result reads 


B. Finite density of grains 


Consider a sample of a rectangular shape with a width, 
L, and thickness, d, (L ^ d). As the injected spin cur¬ 
rent flows through the cross section, the voltage builds 
up between the edges z = TI//2. The easiest way to 
calculate this voltage is to sum the contributions of in¬ 
dividual dipoles. If a grain is located at a point with 
coordinates (Xi,^^,^^), see Fig. then the potential dif¬ 
ference between the edges, created by an induced dipole 
reads 


V{Xi,yi,Zi) = 




[x^ + + (f - Zi)^] 2 


(§■ + 


{x^ + yf + (.^ + z,y]i 


), (19) 


where Vc is given by Eq. (16). In calculating the effective 


inverse spin Hall voltage the summation over dipoles is 
replaced by integration 


* 


{yo)=n dx dy dzV{x,y,z), 

J-i J-yo J-^ 


( 20 ) 


Kf(2/o) = 2ndP, 


In 




— In 



( 21 ) 


At small distances from the injection point, d<^yo<^L, 
the first two terms in Eq. (21) dominate. The second 
logarithm describes a gradual increase of with y^. 
At large distances, y^ the second and the third 

logarithms combine into ln{L/d) leading to the result 


Ktf (oo) = 4n(iPe In (—) 


24e(na^)(iln(^)7 


(^in T 2(Tq^^) 


2Do 


Di, 


Q^i(f) 

An(f) 




( 22 ) 


Note, that for highly conducting grain, both factors in 
denominator do not depend on the characteristics, 
and T^out, of the matrix. In this domain depends 
strongly on the relation between A and a. Overall, Eq. 
(21) describes the growth and subsequent saturation of 


the inverse spin Hall voltage. 


III. MAGNETIC-FIELD DEPENDENCE 


where yo is the distance from the point at which voltage 
is measured to the point of spin-current injection. Natu¬ 
rally, the replacement of the sum by integral is justified 
when nL‘^d ^ 1. The integration over y is straightfor¬ 
ward. Subsequent integral over diverges logarithmically 
at ^ = L/2 and z = —Lj^. This divergence should be 
cut off at ( 2 ( ± I//2) ^ d. Then the integration over x 


The behavior of with position, becomes non¬ 
trivial in the presence of magnetic field directed along 
the ^-axis, a somewhat similar effect was pointed out in 
Ref. |27l If the magnetic field is weak, so that the Lar- 
mour frequency, uJl, is much smaller that and much 
smaller than ^ which is the inverse diffusion time 

through the grain, then the effect of magnetic field on 
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generation of electric dipole can be neglected. Instead, 
the field affects only the polarization in the spin current 
incident on the grain. This allows one to use the result 
Eq. (16) in calculation the cji^-dependence of EJf. 

Outside the grains, the polarization components, Px 
and P^, satisfy the system of equations: Pout^^ + 

^lPz = 0 and - ^lPx = 0. Assuming that 

at the point of injection the polarization was along x, we 
find 


Px{y) =Pa;(o)cos 
Pz{y) =Px{0)sm 


Y ^ 

Yy_ 

exp 


)b_ 


^ ‘^Pout ' 

Y > 

Yy 

exp 

■ / > 

Yv 




(23) 


Suppose that a grain is positioned at ^ = ^o- Then the 
induced dipole moment will be a vector orthogonal to 
polarization with components 


'PziVo) 


PxiVo) 

Px{0) 


Pc, PxiVo) 


PziVo) 


Pc 


where Pc, given by Eq. (16). 


PxiO) 

(24) 

is proportional to the mag- 
ig, which does not change 
To proceed further. 


nitude of the spin current, 
in the presence of magnetic field, 
we notice that only P; 2 -component of the induced dipole 
moment contributes to the buildup of and should be 
substituted into Eq. (19) instead of Pc. We first perform 
integration over z and x. The remaining integral over y 
takes the form 


r 

y!u{yo,^L) = 2ndV^ / 
Jo 


dy 


\/{y-yof Pd? s/{y 


Eor uJl = 0 Eq. ( |25[ ) reproduces the limiting cases of 

L, we 


(26) 


Eq. ( |21[ ). With characteristic distance being - 
conclude that characteristic magnetic field is 

1 Pon. 


^L = 7^ = 


L 2 


which is a natural scale at which the diffusion time 
through a square with a side L is equal to the Larmour 
period. Simple asymptotic expressions for can be 
obtained in the domain uJl ^ uJli when the second term 
in the integrand can be neglected: 

/ \ 1/2 

(i) d yo (Pout/^ lJ • In this limit, the log- 

divergence at large y is cut off at ^ ^ (Pout/^i.)^^^, and 
we get 


Ejf (^0,^1.) = 2ndPcln 
1/2 


(/yi) (27) 


/ \ 

(ii) ^0 ^ [Pout/^L] • We can now neglect y com¬ 
pared to yo in the square brackets. Then the integration 
can be easily performed yielding 


Kf{yo.uj,)=ndV^ 




(28) 


The asymptotes Eq. (27), (28) do not cover the entire 
domain of At the crossover field oJl ~ Doct/yl- Eq. 
(27) exceeds Eq. ( [^ by a large factor ^ ln(^o/c^)- As the 


magnetic-field dependence of voltage is plotted numeri¬ 
cally, see Eig. it appears that in the intermediated do¬ 
main the ISHE voltage exhibits two sign reversals. This 
means that the oscillations in Eq. (23) do not average 


out completely after integration over the positions of the 
spheres. 


1 

yo?PL^ 



exp 



(25) 



FIG. 5. (Color online) The effective spin diffusion length in 
the units (37r^na)“^/^ is plotted from Eq. (36) for ratios 
Pout/Pin: Pout/Pin = 1 (grccu), Pout/Pin = S^urplc) , and 
Pout/Pin = 10 (blue). Note the saturation of Aetf at small A . 


IV. EFFECTIVE SPIN DIFFUSION LENGTH 


There are two reasons why the effective spin-diffusion 
length of the mixture exceeds A. The first reason is obvi¬ 
ous: the grains are sparse and there is no spin relaxation 
in between the grains. The second reason is much more 
subtle and becomes important when A is much smaller 
than the grain radius. Namely, the rate of the spin re¬ 
laxation at the grain surface is suppressed. Formally, 
this suppression, illustrated in Fig. follows from the 
behavior of polarization inside the grain 
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PUr) = - 


3a /*i(t) 






D,^M\h(f). 

_ 3an(x) _ 

2Do„tii(f)+A„K(f) 


i, sinOe^ 


(29) 


It is seen from Fig. that, for A = a/12, the radial 
distribution of Pin(^) not only falls off rapidly from the 
surface towards the center, but its value at the surface is 
small. Physical origin of this smallness is elucidated in 
the Appendix. 

While our goal is to find Agff, in order not to deal with 
boundaries we first calculate the effective spin relaxation 
time of the mixture. Spin relaxation takes place only 
inside the spheres. If at time t = 0 the polarization inside 
the sphere is distributed according to Eq. (29), then the 


rate of decay of this polarization is given by the integral 
over the volume of the sphere 



(30) 


Using the explicit form, ii(x) = 

(xcosh(x) — sinh(x))of the modified spherical 
Bessel function, the integral can be evaluated, and the 
result can be cast in the form 


Substituting Eq. (33) into Eq. (34), we readily find 




'^eff ■ 


37r^na^B^^,Fm' 


(35) 


Note that the product in the numerator is equal to A^. 
We can now use the expression for the effective relaxation 
time to find the effective spin-diffusion length 


'^eff ^ o 


[37r^na^E(^)] 


1/2- 


(36) 


Let us trace the decrease of Aeff as the spin-diffusion 
rate inside the sphere gradually decreases. Eor A ^ a, 
the function F{x) can be replaced by E(0) = 4 (2D^t+^- ) ‘ 
Thus the enhancement of the spin-diffusion length due to 
patterning the SO material into granules is ^ (na^)“^/^. 
In the opposite limit X a we have F{x) ~ 1/x^. This 
leads to the unexpected conclusion that in this limit Aeff 
saturates at the value ^ origin of this sat¬ 

uration is suppression of polarization at the surface, the 
effect discussed above and further elaborated on in Ap¬ 
pendix. 


V. DISCUSSION 


A„t, VA/’ 


(31) 


where the dimensionless function F(a;) is defined as 


F(a;) = 


X sinh(a;) — 2 cosh(a;) + 2 


An 


1) [x cosh(x) — sinh(x)]x -h x^ sinh(x) 

(32) 

The result Eq. ( [M] ) can be also expressed through 
the polarization outside the sphere by replacing by 
PontDout see Eq. One has 


R = 


37r^a^PoutA 


out -'-^out 


D:„T, 




(33) 


(i). The two main results of the present paper are 
Eqs. (21), ([^ and Eq. (36) for the effective inverse 


spin Hall voltage and the effective spin diffusion length 
of the mixture. It is convenient to cast Eq. ( [2^ in the 
form of the relation between the effective spin-Hall angle, 
of the mixture and the spin-Hall angle, of the 
material of the grain. The spin Hall angle is defined as 
the proportionality coefficient between the charge and 
spin current densities, more precisely, ?’ = 0^^(2e/h)js. 
Then Eq. ^ takes the form 


0 


SH 

eff 


12(na3)ln(^) 

2Dout I ’ 

An Ail(f) 


(37) 


In the absence of spin current, the spin relaxation in¬ 
side the spheres causes the time decay of the spin po¬ 
larization in the medium between the spheres. This is 
because diffusing carriers eventually “hit” a sphere. Con¬ 
sider an interval (^o — + ^)5 and assume that there 

are hard walls at the ends, so electrons do not flow in 
or out. Then the initial net polarization, Pout{yo)^y, in¬ 
side the interval will decay with some effective rate 
To find this rate, we substitute the the two-dimensional 
density of spheres in the interval, nSy, into the balance 
equation 


= nSyR. 

Uff 


(34) 


where we assumed ^ Essentially, the propor¬ 

tionality between 0^^ and 0^^ is determined by a “volume 
factor”, na^. Note, however, that 0^^ is the characteris¬ 
tics of a homogeneous S0-film, only as long as the film 
thickness, re, is much smaller than A. Eor re ^ A, 0^^ 
falls off as X/w. At the same time, the decay of 0^^ with 
yo sets in only when yo exceeds the effective spin diffusion 
length of the mixture. This length is much bigger than 
A, as it was shown in Section IV. 

(ii) Note that, strictly speaking, Eq. (36) describes Aeff 
only within a numerical factor. This factor was lost as we 
replaced by PoutDont/a, assuming that the first term 
in Eq. ^ dominates. In fact, precisely at r = a, the two 
terms almost cancel each other. Indeed, substituting the 
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expression for into Eq. ® , we can cast it in the form 


P = — 

nnt 


3a-^ 






\.0e. 


ADoutnCf) 

(38) 

a, the expression in the 
, and thus, is much smaller 
^ a the compensation of 


In the limit X a and r = 
brackets is equal to SAZ^o^t/An 
than a. However, for bigger r 
the first two terms does not take place, and the relation 
is ~ -Pout-Do„t/a holds. 

The suppression of near the surface of the sphere, 
expressed by Eq. ( [3^ , is the reason why Aeff saturates 
when A ^ 0, see Eig. Loosely speaking, strong re¬ 
laxation “repels” the spins from the boundary, which, 
in turn, slows down the effective relaxation. The above 
physics is quite general. To illustrate it, in the Appendix 
we consider a model example of diffusion of particles in 
the presence of an absorbing trap and demonstrate that, 
with increasing the absorption rate, the concentration of 
particles vanishes at the position of the trap. 

(Hi). It is instructive to compare our result Eq. 


( 22 ) 


with the expression for the perturbation of spin cur¬ 
rent flowing in a normal metal around a ferromagnetic 
spher^^. Rather that the SO coupling in our case, the 
difference of spin-up and spin-down carriers in Ref. |23 
is caused by the difference of their conductivities inside 
the ferromagnet. As a result, the induced dipole mo¬ 
ment in our case is normal to the spin current, while the 
induced “spin dipole moment”!^ is along the spin cur¬ 
rent. Other than that, the two expressions resemble each 
other. There is, however, an important difference. If 
the conductivity of the ferromagnetic spher^^ is much 
higher than the conductivity of the surrounding normal 
medium, then the perturbation of the spin current is sup¬ 
pressed (resistance mismatch). On the contrary, for the 
inverse spin Hall effect, the bigger is the ratio An/T^out, 
the stronger is the modification of the spin current out¬ 
side the sphere. 

{iiii). Eor a quantitative example of the effect of granu¬ 
larity on the effective parameters of the mixture, assume 
that the density of the SO granules is na^ = 10“^, while 
the spin diffusion length in the material of the granule is 
A = 0.2a. Compared to the geometry in Eig. [^with no 
granularity we “lose” 100 times in the inverse spin Hall 
voltage. At the same time, we gain in Aeff. Substituting 
A = 0.2a into Eq. 
find Aeff = lOA. 


(36), and assuming ^ 


we 


VI. CONCLUDING REMARKS 

1. Eor experimentally verifying our theoretical results, 
composites of SO and no SO materials can be prepared 
using a variety of widely available fabrication techniques. 
Eor example, in Ref. [29l authors used a pulsed laser de¬ 
position technique to prepare a composite comprising of 
gold nanoclusters embedded in ZnO matrix. In Ref. [30l 


a self-assembly approach was used to fabricate a compos¬ 
ite comprising of nickel nanoclusters embedded in amor¬ 
phous AI 2 O 3 matix. In Ref. [3TJ a nanofabrication ap¬ 
proach was employed to prepare a magnonic crystal com¬ 
prising of cobalt nanodots embedded in a permalloy film. 
Similar approaches can be used to prepare the desired 
composite structures of SO and no SO materials, say Pt 
or An nanodots (with large 0^^) embedded in films of 
low 0^^ materials (such as copper, molybdenum or even 
semiconductors like silicon). 

2. Erom device point-of-view, an obvious way to en¬ 
hance the spin diffusion length would be by creating a ID 
structure of alternating SO and no SO layers To achieve 
this, however, the thickness of the SO-layer should be 
smaller than A. Conversely, in granular system the en¬ 
hancement takes place when X a. This is because the 
spin current can flow around the spheres. 

3. In numerous spin-pumping experiments, see e.g. 
Refs. [ 4 IIT 7 I the measured quantity, is proportional 
either to ^^^A, when the thickness of non-magnetic ma¬ 
terial is much bigger than A, or simply to 0^^ in the op¬ 
posite limit. A comprehensive list of experimental values 
of 0^^ and A for a number of heavy metals can be found 
in Ref. [ 32 I This list indicates that, while, separately^ 
and A vary within wide ranges, the range of change of 
their product is much narrower, see also Ref. [33| 

Overall, there is still experimental ambiguity in ex¬ 
tracting the intrinsic SO parameters of materials from 
the experiment. In this regard, granularity can offer a 
help, by bringing a new spatial scale, the radius of the 
grain, a. As shown in Eig. the value Aeff depends very 
strongly on the relation between A and a. 

4. In a specific case of a semiconductor ZnO the 

inverse spin Hall effect was studied both in pumping 
experimenfP^ and directly by measuring the nonlocal 
voltag^^. In both measurements the value 0^^ w as found 
to be anomalously big , compared e.g. to . It has 

recently been showiP^ that the value 6^^ in ZnO can be 
tuned very sensitively by changing the oxygen ambient 
under which it is growrP^. Eilms prepared under high 
oxygen rich environment showed a large value for 0^^ 
(^ 0.1), while the films prepared under low oxygen am¬ 
bient showed an order of magnitude lower value of 0^^. 
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Appendix A 

Consider a diffusion in one dimension. If at time t = 0 
the distribution of particles is a (5-peak, i.e. n(x, 0) = 











n{x,to)[arb.units] 


n{x, to) [arb.units] 



FIG. 6. Diffusive spreading of the initial particle distribution 
n(x,0) = S{x) in the presence o f an a bsorbing trap located 
at a: = 0 is described by Eq. (AlO). Shown is n{x,t) at 
a fixe d time, to, for different absorption efficiencies, k, Eq. 
(A13). The more absorbing the trap is, the deeper is the dip 


at the origin, and the slower is the decay of t he net number 
of particles at long times, as follows from Eq. (A12). 


S{x)^ then it spreads with time as 

where D is the diffusion coefficient. 

Suppose now, that an absorbing trap is placed at the 
coordinate origin. Then the spreading is governed by the 
equation 


dn ^d‘^n a6(x)n 

— = -D —^ H- 

dt dx^ r. 


(A2) 


where a is the size of the trap, and is the absorption 
rate. Then the time-dependent concentration, n( x,t) , 
can be expressed through the Green function of Eq. ( |A2| ) 

n{x^t) = J dx'G{x^x'^t)n{x\0) (A3) 

It is convenient to present the Green function in terms 
of eigenfunctions, of Eq. (A2), which satisfy the 
Schodinger-like equation 


, aS{x) 2 , 

- H-= k 

OX^ Ts 

Then the expression for G(x,x’,t) reads 

G{x,x',t) = y^ tpk{x)‘>pk{x') ex.Tpj-Dk'^t). 


(A4) 


(A5) 


The second term in Eq. ( |A4[ ) plays the role of delta- 
1 car 

d'lpk 


potential barrier, and causes the discontinuity of the 
derivative of tpk 


d'lpk 

dx 


a:= 0 + 


dx 


^ (A6) 

x=o- JJr^ 


The normalized solutions, which satisfy Eq. (A4) 

have the form 


1 


'ipk{x) = cos(A:|a:| + ipk), 


where the phase (pk is found from the condition 

a 

idxipk = ~ 


2DtX 


(A7) 


(A8) 


imposed by Eq. (A6). 


Upon substituting Eq. (A5) into Eq. 


and using 


the initial condition, we arrive at the final result 


n(x,t) = — 


TT Jo 


cosk\x\ + 
dk- 




^ + ^2DT.k) 


(A9) 

It is now convenient to introduce a dimensionless coor¬ 
dinate X = and the dimensionless time de pen- 

dent parameter I = In new variables Eq. (A9) 

assumes the form 


n(x, t) = 


27t Dr. Vt 


7 cos \x 
= / dq - 

i Jo 


\q^ Vtqsm\x\q _ 2 




(AlO) 

We see that the characteristics of the trap, a and r^, 
enter only into rescaling of time. In Fig. [^we plot Eq. 
(AlO) for four different i. It is seen from Fig. [^that. 


with time, the density n(0, t) at the origin develops a 
dip. The smaller is r^, i.e. the more absorbing is the 
trap, the sharper is the dip. This conclusion also follows 
from the long-time asymptote of n(x, t), when we neglect 
q^ in denominator compared to t. Then the integration 
yields 


n(x, t) 


t>i 


S^/7^ Dr si 

/ a\x 


\2Dt, 


-h l^e "4 


(All) 


This asymptote indicates that the ratio of concentrations 
at half-width, x ~ (I)t)^/^, and at the origin is ^ 
i.e. the dip is deep. 

The next question we ask ourselves is how the to¬ 
tal number of particles f dx n(x, t) decreases with time. 
Upon integration Eq. (AlO), we get 


/ oo 

dx n(x, t) 

-00 

=Erfc(f2) exp(t) = < 


i-ft(A)’. ‘«i. 


*»i. 


(A12) 


where Erfc(s) is the complementary error function. It is 
seen from Eq. ( |A12 ) that the change of the decay rate 
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dN/dt takes place at t ^ 1. This change is caused by the 
development of the dip. Indeed, for t ^ 1 the decay rate 
falls off with time as 

Overall, we conclude that the spreading of the parti¬ 
cle density in the presence of a trap is governed by a 
dimensionless parameter 


then, n{x^t) will develop a dip at x = 0 and the decay of 
the net number of particles will proceed even slower. For 
large efficiency, the dip will developed early, namely at 
t ^ <C Ts, after which time the decay of N(t) will 

be governed by the value of n(0, t) of the concentration 
at the dip. 


a 


(A13) 


which is the dimensionless efficiency of absorption by the 
trap. If this efficiency is small, the spreading will proceed 
as in the absence of the trap for most of the time, until 
the concentration at x = 0 becomes really small. Only 


Formation of a dip in our model problem puts into 
a general perspective the behavior of the effective spin 
diffusion length in the system of the SO grains. In the 
limit A ^ a, see Eq. (36), the value Agff saturates because 
the polarization near the boundary gets suppressed as a 
result of the development of a local minimum. 
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